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The compressibility of nuclear matter has received significant attention in the last decade and a variety of 
approaches have been employed to extract this fundamental property of matter. Recently, significant differences 
have emerged between the results of relativistic and non-relativistic calculations of breathing mode giant monopole 
resonance (GMR). This is due to a lack of understanding of the dynamics of GMR and of its exact relationship 
to the compression modulus of the infinite nuclear matter. Here, I present an alternative approach based upon 
nuclear shell effects. The shell effects are known to manifest experimentally in terms of particle-separation energies 
with an exceedingly high precision. Within the framework of the non-relativistic density-dependent Skyrme theory, 
it will be shown that the compressibility of nuclear matter has a significant influence on shell effects in nuclei. I 
will also show that 2-neutron separation energies and hence the empirical shell effects can be used to constrain 
the compressibility of nuclear matter. 



Introduction 

Shell effects manifest themselves in terms of the existence 
of magic numbers in nuclei. Inclusion of spin-orbit cou- 
pling in the shell model is known to account for the magic 
numbers 1 -'. Shell effects are also found to appear more 
subtly e.g. as anomalous kink in the isotope shifts of Pb 
nuclei. It was shown that in the relativistic mean-field 
(RMF) theory, this kink appears naturally due to the in- 
herent Dirac-Lorentz structure of nucleons 2 ). By includ- 
ing an appropriate isospin dependence of the spin-orbit 
potential in the density dependent Skyrme approach, the 
anomalous kink in the isotope shifts in Pb nuclei could 
be produced 3 ^ . 

The shell gaps at the magic numbers are known to pro- 
■ duce important effects. Consequences due to the shell 
gaps are broadly termed as shell effects. The clearest 
manifestation of the shell effects can be seen in the pro- 
nounced dip in neutron and proton separation energies at 
the magic numbers in figures by Wapstra and Audi 4 ) on 
the experimental separation energies all over the periodic 
table. Recently, shell effects have become focus of much 
attention 5 ' 6 ^. This is due to the reason that the shell 
effects in nuclei near the drip lines constitute an impor- 
tant ingredient to understand r-process abundances and 
heavy nucleosynthesis 7,8 ) . 

In this talk, I focus upon 2-neutron separation energies 
to demonstrate that the shell efffects show a significant 
dependence upon the compressibility of nuclear matter. 
Using empirical data on the 2-neutron separation ener- 
gies and thus implicitly the empirical shell effects, I will 
show that these data can be used to constrain the com- 
pressibility of nuclear matter. As the focus of this confer- 
ence shows, the compression modulus K is an important 
fundamental property of the nuclear matter. It repre- 
sents a cardinal point on the behaviour of equation of 



state (EOS) of the nuclear matter. 

The compressibility of the nuclear matter has received 
a significant attention in the last decade and various 
approaches 9 " 11 ) have been employed to extract the com- 
pression modulus. However, an unambiguous determina- 
tion of the compression modulus of the nuclear matter 
has remained difficult due to a lesser sensitivity of the gi- 
ant monopole resonance data to the compressibility 12,13 ). 



The Density-Dependent Skyrme Theory 

Here I employ the density-dependent Skyrme theory 14 ) 
to examine the shell effects. The Skyrme approach has 
generally been found to be very successful in provid- 
ing ground-state properties of nuclei. The energy den- 
sity functional used in the Skyrme approach is the stan- 
dard one as given in Ref. 15 ). The corresponding energy 
per nucleon for the symmetric nuclear matter (with the 
Coulomb force switched off) is given by 



{E/A)oo = kpV\l + Pp) + h p + ^t 3 p 1+a 

8 16 



(1) 



where k — 75 MeV.fm 2 . The constant (3 is given in terms 
of the constants t\ and *2 and x 2 of the Skyrme force by 



h 2 4 
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(2) 



The constant (3 is also related to the effective mass m* 
and the saturation density po by 



m/m* = 1 + f3po 



(3) 
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The incompressibility (or the compression modulus) of 
the infinite nuclear matter is given as the curvature of 
the EOS curve and can be written as 



K 



9Po 
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dp 2 



(4) 
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The parameters to, t^ and a are usually obtained from 
the nuclear matter properties. The other parameters t\, 
i 2 and various x parameters are obtained from fits to 
properties of finite nuclei. The strength Wo responsible 
for the spin-orbit interaction is obtained by reproducing 
the spin-orbit splittings in nuclei such as 16 and 40 Ca. 



I have fixed x\ and X2 at zero for convenience. Given 
a value of (3, the coeffients t\ and t 2 in cq. (2) are ob- 
tained from the fits to the ground-state binding energies 
of finite nuclei. Thus, I have obtained various Skyrme 
parameter sets with incompressibility values K = 200, 
220, 249, 270, 305, 327, 360 and 393 MeV, respectively. 
These forces encompass a broad range of physically plau- 
sible values of K. The nuclear matter properties of these 
forces are shown in Table 1. It can be seen that the 
force with K = 270 MeV was obtained for the sake of 
interpolation. Table 2 shows the total binding energies 
of some key nuclei obtained with the HF+BCS approach 
using the various Skyrme forces. It can be seen that all 
of these forces reproduce the total binding energies of the 
key nuclei from 16 to 208 Pb very well. 



The Skyrme Forces 

In order to show the effect of the incompressibility of 
nuclear matter on shell effects and consequently also on 
ground-state properties of nuclei, I have constructed a 
series of zero-range Skyrme forces. Experimental data 
on ground-state binding energies and charge radii of key 
nuclei such as 16 0, 40 Ca, 90 Zr, 116 Sn, 124 Sn and 208 Pb 
are taken into account in the least-square minimization. 
The ground-state properties in the Skyrme theory are 
calculated using the Hartree-Fock method. 

With a view to vary K over a large range, the correlation 
between K and the saturation density is implicitly taken 
into account. This correlation has been summarized for 
Skyrme type of forces in Fig. 4 of Blaizot 9 ' . Accordingly, 
there exists an inverse correlation between the saturation 
density (Fermi momentum) and the compression modu- 
lus. Keeping this in mind, I have varied the saturation 
density po over the range 0.140 - 170/m- 3 in steps of 
.005/to~ 3 as shown in Table 1. 

The effective mass is fixed at 0.79. This value is required 
for being consistent with giant quadrupole resonance en- 
ergies of heavy nuclei 15 ). Eq. (3) then provides the value 
of (3 which can be used in eq. (2) to connect the coeffi- 
cients t\ , ti and x 2 ■ In order to do a systematic variation 
of the nuclear compressibility, I have kept the saturation 
binding energy fixed at a value of —16.0 MeV. This value 
is consistent with most nuclear mass models 16 ) and is 
close to physically acceptable values. 

The coefficients a, to and £3 are determined by the eqs. 
(1), (4) and the saturation condition. However, we have 
allowed to and £3 to vary so that a fit to finite nuclei 
provides an incompressibility that is consistent with the 
inverse correlation. On the other hand, in an exhaustive 
computational exercise it is found that the violation of 
the above correlation results in bad fits for ground-state 
binding energies of key nuclei. 



Two-neutron Separation Energies and the Compressibility 

I have selected the chain of Ni isotopes in order to probe 
the shell effects. As most Ni isotopes are known to be 
spherical and the experimental binding energies are known 
over a large range of Ni isotopes, the chain of Ni isotopes 
serves as an ideal test bench to probe the shell effects. 
Even-mass Ni isotopes from A=52 to A=70 are consid- 
ered. This includes the neutron magic number N=28 
(A=56) where we intend to investigate the shell effects 
due to the major shell closure. 

With a view to obtain the ground-state properties of nu- 
clei, we have performed spherical Skyrme Hartree-Fock 
calculations in coordinate space. Herein pairing is in- 
cluded using the BCS formalism with constant pairing 
gaps. Since nuclei under focus such as 56 Ni and 58 Ni are 
not far away from the stability line, the BCS scheme suf- 
fices to be a suitable mechanism for the pairing. Ground- 
state binding energies and rms charge and neutron radii 
for the Ni isotopes are calculated using the various Skyrme 
forces. The binding energies of nuclei are used to obtain 
the 2-ncutron separation energies as 

S 2n (Z,N) = B(Z,N)-B(Z,N-2), (5) 

where B represents the total binding energy of a nucleus. 

As mentioned earlier, effects due to the shell gaps affect 
various nuclear properties besides the particle separation 
energies. First, the calculated charge radii of Ni isotopes 
are shown in Fig. 1. The charge radii show an increase 
with the incompressibility. Evidently, forces with the in- 
compressibility about 200-250 MeV show a little disper- 
sion in the values of the charge radii. However, above 
K ~ 300 MeV, there is a clear increase in the charge 
radii as a function of K and forces with a large incom- 
pressibility yield large values of charge radii for any given 
isotope. Inevitably, a large value of K hinders synthe- 
sis of nuclei with a radial extension which is diminutcd 
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(compressed) contrary to that with a lower values of K. 
The empirical rms charge radii of 58 Ni, 60 Ni, 62 Ni and 
64 Ni are taken from the compilation of Ref. 17 ) by solid 
circles. The values are 3.776, 3.815, 3.846 and 3.868 fm, 
respectively, obtained by folding the proton density dis- 
tributions with the finite size 0.80 fm of protons. These 
results derive from experiments on muonic atoms. The 
rms charge radii deduced from the precision measure- 
ments in muonic atoms (taken from the recent compi- 
lation of Ref. 18 ) are 3.776, 3.813, 3.842 and 3.860 fm, 
respectively, which are almost identical to the previous 
compilation 17 ) . 

The empirical charge radii shown in Fig. 1 encompass the 
theoretical curves between K=270-327 MeV. Whereas 
the charge radius for 58 Ni points towards K=270 MeV, 
that for the heavier Ni isotopes is in the vicinity of K=327 
MeV. Since the charge radii deduced from various meth- 
ods have been obtained with a significant precision, these 
data could, in principle, be used to constrain the incom- 
pressibility. However, as there is still some model depen- 
dence in extraction of the charge radii, I would like to 
withhold any conclusions from this figure. On the other 
hand, the general trend of the experimental data points 
to a higher value of the compression modulus. 

Figure 2 shows the corresponding rms neutron radii. The 
neutron radii show a monotonous increase with the mass 
number, with the exception of a slight kink at the magic 
number N=28. The change in the r n values with K is 
similar to that for the charge radii, i.e. for lower val- 
ues of K, there is a very little change in r n with K. A 
significant change in the neutron radii, however, can be 
seen with large K values. The empirical rms neutron 
radii obtained from 800-McV polarized-proton scatter- 
ing experiment 17,19 ^ for the isotopes 58 Ni and 64 Ni are 
shown in the figure by the solid points. The values lie 
between the curves for K=300 and K=327 MeV. 

The 2-neutron separation energy S 2 n for the Ni isotopes 
obtained with the various Skyrme forces are shown in 
Fig. 3. For the sake of clarity of the presentation, I have 
selected a set of the Skyrme forces in the figure. The 
dramatic fall in the 52„ curves is seen conspicuously for 
the nucleus just above the N=28 magic number. Such a 
kink in the S^n values signifies the presence of the shell 
effects which arise from closure of a major shell. All the 
Skyrme forces produce such a kink in Fig. 3. However, 
the slope of the kink between A=56 (N=28) and A=58 
(N=30) changes from one force to the other (as a func- 
tion of the incompressibility K). The difference in the 
S2n values of 56 Ni and 58 Ni can be taken as a measure of 
the shell effects. The force with the lowest incompress- 
ibility (K = 200 MeV) gives this difference as 6.59 MeV 
with a slope which is minimum amongst all the curves. 
This implies that the shell effects due to this force are 
the weakest one. As the K value increases, the corre- 
sponding difference shows a smooth increase and the en- 
suing steepness increases gradually. For the force with 



the highest incompressibility (K = 393 MeV), the dif- 
ference amounts to 9.46 MeV. This difference is higher 
than the experimental difference 20 - 1 of 8.37 MeV. It im- 
plies that the shell effects show a strong dependence on 
the compression modulus K. Notwithstanding this cor- 
relation, the difference in the experimental values of Si n 
can serve as a calibration for K. Consequently, the exper- 
imental difference and the slope is closest to the curve for 
K=270 MeV. Thus, on the basis of the empirical Sm val- 
ues, I infer that the incompressibility of nuclear matter 
lies in the neighbourhood of K ~ 270 MeV. 



Shell Effects at the Neutron Drip-Line 

Shell effects near the drip lines have been a matter of dis- 
cussion of late 5 ' 6 ). It has been shown earlier that within 
the non-relativistic approaches of the Skyrme type, the 
shell effects near the drip lines and in particular near 
the neutron drip line are quenched 5 ). This is contrary 
to that shown in the RMF theory by Sharma et al. 6 ), 
where the shell effects were observed to remain strong. 
Before such a debate is resolved, it is important to set- 
tle the issue of the shell effects near the stability line. 
Fig. 4 shows the 2-neutron separation energies for the Ni 
isotopes using the force SkP within the Hartree-Fock Bo- 
golieubov approach 21 ). The force SkP has been found to 
be successful in reproducing the ground-state properties 
of nuclei. The comparison of the SkP results with the ex- 
perimental S2n values shows that empirically there exist 
strong shell effects than predicted by the Skyrme force 
SkP. The compression modulus of the force SkP is about 
200 MeV. Thus, the behaviour of the force SkP about 
the weak shell effects is consistent with the results from 
forces with low value of the incompressibility in Fig. 3. 
This suggests again that the shell effects in Ni are com- 
mensurate to a higher value of the compression modulus 
of the nuclear matter as inferred above. 

How the shell effects behave along the neutron drip line 
as a function of the compression modulus can be visu- 
alized in Fig. 5. Here I have chosen the chain of Zr 
isotopes. This chain includes nuclei on both the sides 
of the magic number N=82 (A=122). The total bind- 
ing energies obtained in the Hartree-Fock calculations 
with the various Skyrme forces are shown. For nuclei 
from A=116 (N=76) to A=122 (N=82), the total bind- 
ing energy shows a steep increase in the value with all the 
forces. The forces with lower K show a stronger binding 
in general as compared to those with higher K, which is 
as expected. For nuclei above A=122, there is a strik- 
ing difference in the way the binding energies progress 
with mass number. For the force with K = 200 MeV, 
the binding energy of nuclei heavier than 122 Zr shows an 
increase in the value with the mass number, implying 
that a further sequential addition of a pair of neutrons 
to the N=82 core does contribute to the binding energy. 
As the K value increases, the binding energy contribu- 
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tion from neutrons above N=82 starts diminishing. For 
K values 300 MeV and above, a stagnation is observed 
in the binding energies. Such a behaviour implies that 
the shell effects become stronger as the incompressibil- 
ity increases. This is again an indication that the shell 
effects are strongly correlated to the incompressibility of 
the nuclear matter. Consequently, the incompressibility 
inferred from the ground-state data along the stability 
line suggests that the shell effects near the neutron drip 
line are strong. This is consistent with the strong shell 
effects about the neutron drip line as concluded in the 
RMF theory 6 ). 

For nuclei near the drip lines, the HFB approach is con- 
sidered to be more suitable than the HF+BCS one. How- 
ever, the effects of the continuum are expected to be 
significant more for nuclei very near the drip line and 
in particular the observables which are most affected by 
coupling to the continuum are the rms radii. The to- 
tal binding energies are known to show a little difference 
in the HF+BCS and HFB approaches. The main char- 
acter of the shell effects is not altered by the inclusion 
of the continuum, as has been pointed out in Ref. 8 -* us- 
ing the Skyrme force SIII. However, it will be interesting 
to see how a judicious choice of a force compatible with 
the experimental data as discussed above, will affect the 
r-process nucleosynthesis. 



Summary and Conclusions 

In the present work, I have shown that there exists a 
correlation between the shell effects and the compress- 
ibility of nuclear matter. The latter is shown to influ- 
ence the shell effects significantly. Consequently, the 2- 
neutron separation energies show a strong dependence 
on the compression modulus. The ensuing correlation 
then provides a calibration for the compression modulus 
of nuclear matter. It is shown that using this correlation, 
the empirical data on the 2-neutron separation energies 
and thus implicitly the shell effects can be used to con- 
strain the compressibility of nuclear matter. This proce- 
dure allows me to conclude that the incompressibility K 
should be within 270-300 MeV. This conclusion is also 
consistent with the experimental data on charge radii. 
A natural consequence of the present study is that the 
shell effects near the neutron drip line are predicted to be 
strong. This is in conformity with the strong shell effects 
observed at the neutron drip line in the RMF theory. 
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Tabic 1. The nuclear matter properties of the various 
Skyrme interactions. The saturation energy E/A has 
been fixed at —16.0 MeV and the effective mass m* 
has been kept at 0.79. 



Force 


Po {fm a ) 


K (MeV) 


a sym (MeV) 


I 


0.140 


393 


38.0 


II 


0.145 


360 


32.4 


III 


0.150 


327 


32.5 


IV 


0.155 


305 


32.3 


V 


0.160 


249 


32.5 


VI 


0.165 


220 


32.0 


VII 


0.170 


200 


31.8 



Table 2. The total binding energy (MeV) of nuclei 
obtained from the HF+BCS calculations with the 
Skyrme interactions. The empirical values (exp.) are 
shown for comparison. 



K (MeV) 



Nucleus 


exp. 


393 


360 


327 


305 


249 


220 


200 




-127.6 


-127.5 


-127.5 


-125.5 


-127.5 


-127.5 


-127.6 


-127.7 


40 Ca 


-342.0 


-342.5 


-342.4 


-342.2 


-342.0 


-342.4 


-342.5 


-342.4 


90 Zr 


-783.9 


-786.6 


-786.0 


-785.6 


-785.2 


-785.3 


-784.6 


-784.3 


116 Sn 


-988.7 


-989.2 


-988.1 


-987.0 


-986.1 


-986.1 


-985.5 


-985.4 


120 Sn 


-1020.5 


-1021.3 


-1020.1 


-1019.6 


-1019.3 


-1019.1 


-1018.9 


-1018.9 


124 Sn 


-1050.0 


-1051.0 


-1049.7 


-1050.0 


-1050.5 


-1049.9 


-1050.1 


-1050.3 


208p b 


-1636.5 


-1636.5 


-1635.7 


-1636.0 


-1636.5 


-1636.2 


-1637.2 


-1637.0 



Figure Captions 

Fig. 1. Charge radii for Ni isotopes with various Skyrme 
forces. The experimental values for a few nuclei 
from Refs. 17 ' 18 -* are shown by solid dots. 

Fig. 2. Neutron radii for Ni isotopes with various Skyrme 
forces. The experimental values for 58 Ni and 64 Ni 
from Refs. 19 ) are shown by solid dots. 

Fig. 3. Two-neutron separation energy S2 n for Ni iso- 
topes obtained with various Skyrme forces. Com- 
parison with the experimental data from Ref. 20 ) is 
also made. 

Fig. 4. Total binding energy of Zr isotopes near the neu- 
tron drip line calculated with the Skyrme forces. 
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